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ABSTRACT 

The infrared flux method (IRFM) has been applied to a sample of 135 dwarf and 36 giant stars 
covering the following regions of the atmospheric parameters space: 1) the metal-rich ([Fe/H] > 0) end 
(consisting mostly of planet-hosting stars), 2) the cool ( T e g < 5000 K) metal-poor (—1 < [Fe/H] < —3) 
dwarf region, and 3) the very metal-poor ([Fe/H] < —2.5) end. These stars were especially selected to 
cover gaps in previous works on T e g vs. color relations, particularly the IRFM T e g scale of A. Alonso 
and collaborators. Our IRFM implementation was largely based on the Alonso et al. study (absolute 
infrared flux calibration, bolometric flux calibration, etc.) with the aim of extending the ranges of 
applicability of their T e g vs. color calibrations. In addition, in order to improve the internal accuracy 
of the IRFM T e g scale, we recomputed the temperatures of almost all stars from the Alonso et al. work 
using updated input data. The updated temperatures do not significantly differ from the original ones, 
with few exceptions, leaving the T e g scale of Alonso et al. mostly unchanged. Including the stars with 
updated temperatures, a large sample of 580 dwarf and 470 giant stars (in the field and in clusters), 
which cover the ranges: 3600 K< T e g < 8000 K, —4.0 < [Fe/H] < +0.5, have T e g homogeneously 
determined with the IRFM. The mean uncertainty of the temperatures derived is 75 K for dwarfs and 
60 K for giants, which is about 1.3% at solar temperature and 4500 K, respectively. It is shown that 
the IRFM temperatures are reliable in an absolute scale given the consistency of the angular diameters 
resulting from the IRFM with those measured by long-baseline interferometry, lunar occultation and 
transit observations. Using the measured angular diameters and bolometric fluxes, a comparison 
is made between IRFM and direct temperatures, which shows excellent agreement, with the mean 
difference being less than 10 K for giants and about 20 K for dwarf stars (the IRFM temperatures 
being larger in both cases). This result was obtained for giants in the ranges: 3800 K< T e g < 5000 K, 

—0.7 < [Fe/H] < 0.2; and dwarfs in the ranges: 4000 K< T e g < 6500 K, —0.55 < [Fe/H] < 0.25; 
and thus the zero point of the IRFM T e g scale is essentially the absolute one (that derived from 
angular diameters and bolometric fluxes) within these limits. The influence of the bolometric flux 
calibration adopted is explored and it is shown that its effect on the T e g scale, although systematic, 
is conservatively no larger than 50 K. Finally, a comparison with temperatures derived with other 
techniques is made. Agreement is found with the temperatures from Balrner line-profile fitting and 
the surface-brightness technique. The temperatures derived from the spectroscopic equilibrium of Fe I 
lines are differentially consistent with the IRFM but a systematic difference of about 100 K and 65 K 
(the IRFM temperatures being lower) is observed in the metal-rich dwarf and metal-poor giant T e g 
scales, respectively. 

Subject headings: stars: atmospheres — stars: fundamental parameters — infrared: stars 


1. INTRODUCTION 

The effective temperature (T e g) is one of the funda¬ 
mental properties of a star. In principle, it may be de¬ 
termined from the stellar angular diameter 6 and the 
bolometric flux fg Q \ according to: 

/boi = ^T e 4 ff , (1) 

where cr is the Stefan-Boltzmann constant. In practice, 
however, the difficulties and limitations that arise when 
measuring angular diameters restrict the direct deter¬ 
mination of temperatures to a relatively low number of 

1 Affiliated with the Seminario Permanente de Astronomi'a y 
Ciencias Espaciales of the Universidad Nacional Mayor de San Mar¬ 
cos, Peru 


stars. The problem is worse for main sequence stars al¬ 
though interferometric (see e.g. Kervella et al. 2004) 
and transit (Brown et al. 2001) observations have been 
recently performed. 

Even though there are alternative ways to determine 
effective temperatures, all of them require the introduc¬ 
tion of models, which makes the results somewhat de¬ 
pendent on the models adopted. In the range of spec¬ 
tral types F0-K5, the effective temperature may be de¬ 
termined from the spectroscopic conditions of excita¬ 
tion and ionization equilibrium of Fe lines (e.g. Takeda 
et al. 2002a, hereafter T02a), synthetic photometry (e.g. 
Bessell et al. 1998, hereafter B98; Houdashelt et al. 
2000), fitting of the Balrner line-profiles (e.g. Fuhrmann 
1998, hereafter F98; Barklem et al. 2002), and line- 
depth ratios (e.g. Gray 1994; Kovtyukh et al. 2003, 
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hereafter K03). Our approach is the use of the infrared 
flux method (IRFM), which will be described in 1)2 and 
compared to the other methods in ID 

In purely spectroscopic methods, the temperature is 
determined by constraining its value in such a way that 
the abundance of Fe I does not depend on the excitation 
potential of the lower level of the lines used and that the 
average abundance of Fe derived from Fe I and Fe II lines 
are equal (see e.g. T02a). This actually sets a constraint 
on both T e g and logg (surface gravity) so degenerate 
solutions may be found. The technique has been ap¬ 
plied to a wide variety of samples: stars in the thin and 
thick disks of the Galaxy (Bensby et al. 2003, hereafter 
B03), planet-hosting stars (Santos et al. 2004, hereafter 
SIM04), globular cluster stars (Sneden et al. 2004), and 
so on. Evidence is being collected, however, towards the 
fact that the spectroscopic method has systematic errors 
caused by model uncertainties, such as non-LTE effects 
(note also that most spectroscopic studies use 1D-LTE 
classical model atmospheres). For instance, the spectro¬ 
scopic temperatures are ~ 100 K higher than those deter¬ 
mined with the IRFM for planet-hosting stars (Ramirez 
& Melendez 2004b, hereafter RM04b; SIM04), surface 
gravities are ~ 0.1 dex higher than those derived from 
Hipparcos parallaxes for [Fe/H] > —0.7 (Allende Prieto 
et al. 1999, SIM04), and ionization equilibrium of Fe is 
not reproduced by classical model atmospheres neither 
in the Hyades cool dwarfs (Yong et al. 2004) nor in the 
field metal-rich K dwarfs (Allende Prieto et al. 2004, 
hereafter AP04). 

Synthetic spectra and filter transmission functions 
have been used to generate grids of colors in the T e g , 
[Fe/H], log g space using both Kurucz (1979) ATLAS 
models (B98) and Gustafsson et al. (1975) MARCS mod¬ 
els (Houdashelt et al. 2000). Synthetic colors are often 
put into the observational systems with samples of field 
stars (usually with metallicities near solar) with reliable 
temperatures or calibrated color-temperature relations. 
Zero point corrections are finally applied until the T e g 
scale resembles the absolute one. 2 Since this is made 
with field stars of solar metallicity, the zero point cor¬ 
rections are strictly valid only for synthetic colors with 
[Fe/H] ~ 0, and its validity for other metallicities re¬ 
lies on the capability of models to reproduce real stellar 
spectra. Although observations and temperatures de¬ 
rived with other methods are involved in the procedure, 
the synthetic T e g scales retain the basic properties of the 
model atmospheres adopted and are thus useful to test 
them. 

The strength of the hydrogen Balmer lines is very sen¬ 
sitive to the effective temperature below 8000 K and it is 
only weakly affected by the surface gravity. Profile fitting 
is the best approach to an estimate of T e ff with Balmer 
lines, normally within 110 K of random error. Barklem 
et al. (2002) found 5730 ± 70 K for the Sun with up-to- 
date line opacities and a reasonable agreement between 
their temperatures and those obtained with the IRFM for 
T e g > 5800 K. Castelli et al. (1997), however, show that 
the results from Balmer line-profile fitting are strongly 
dependent on the models and the details of the convec- 

2 We refer to the T e fj scale obtained with direct temperatures 
(i.e. stellar angular diameter and bolometric flux measurements) 
as the absolute T e g scale. 


tion treatment. To satisfactorily reproduce the profiles 
in the solar spectrum using ATLAS9 models with over¬ 
shoot, for example, the temperature of the Sun has to be 
increased up to 5900 K. 

With the line-depth ratio technique (Gray & Johans¬ 
son 1991, Gray 1994) very small internal errors (~ 10 K) 
for individual stars may be achieved if several calibrated 
line-depth ratios are used. Even though the temperature 
scale derived with this technique may be less accurate 
due to systematic errors, the line-depth ratios may be 
used to look for very small time dependent T e g variations, 
like those produced by photospheric inhomogeneities and 
starspots. The solar temperature with this method re¬ 
sults to be 44 K below the standard value of 5777 K and 
an empirical correction of 44 K has to be applied (K03). 
Since nothing guarantees that this correction is applica¬ 
ble to the whole range of atmospheric parameters, the 
systematics, if present, are not removed. 

The aim of this work is to extend and improve 
the internal accuracy of the best available empiri¬ 
cal temperature-color calibrations based on the IRFM, 
namely, those given by Alonso et al. (1996b, here¬ 
after AAM96b; 1999b, hereafter AAM99b) 3 for the 
( UBVRI)j(JHK)tcs and Stromgren uvby-/3 systems. 
Melendez & Ramirez (2003) and Ramirez and Melendez 
(2004a, hereafter RM04a) extensions to the Vilnius, 
Geneva, RI(C) an d DDO systems will also be revised. 
In this paper we describe our IRFM implementation 
and test our results. The calibration of the metallicity- 
dependent T e g vs. color relations is the subject of a com¬ 
panion paper (hereafter referred to as Part II). 

This paper is distributed as follows: in 12 the IRFM 
and its present implementation are explained, the sam¬ 
ple and input data adopted as well as their effects on the 
temperatures to be derived are described in [|3J The cal¬ 
culation of the final temperatures and angular diameters 
is given in 0 while the fundamental tests (e.g. compar¬ 
ison with measured angular diameters) and comparison 
with other methods are presented in Conclusions are 
summarized in m 

2. THE IRFM AND ITS PRESENT IMPLEMENTATION 

The IRFM was first introduced by Blackwell & Shallis 
(1977) as a mean to simultaneously determine the effec¬ 
tive temperature and angular diameter of a star. The ba¬ 
sic idea of the method is to use the emergent monochro¬ 
matic flux predicted by models (V’a) and, similarly to 
Eq. Q, given by 

h = (^ff, [Fe/H], log 5 ) , (2) 

where f\ is the corresponding flux as measured on Earth, 
to solve the system of equations and m for the un¬ 
knowns T e ff and 6. The IRFM thus requires /bob [Fe/H], 
log g and f\ to be known a priori. The use of the infrared 
is due to the relatively low dependence of on T e g in 
this region of the spectrum (see e.g. Fig. 1 in Blackwell 
et al. 1979). 

A more practical approach to the IRFM was given 
by Blackwell et al. (1980), who defined the i?-factors 
as the quotients between the total and monochromatic 

3 Hereafter, the series of papers by A. Alonso and collaborators 
are referred to as A AM. 
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fluxes. The comparison of the observational i?-factor 
{Robs = fbol/fx) with a grid of theoretical A-factors 
(i?theo = aT^ s /ip\, note that Atheo is a grid in the T eff , 
[Fe/H], log g space at a given wavelength) allows to easily 
automatize the method. 

The IRFM is often applied using more than one in¬ 
frared wavelength (A), and a different T e g is obtained for 
each A. A critical ingredient of the IRFM is the abso¬ 
lute flux calibration of the standard star from which the 
f\ fluxes of the problem stars are calculated. If /^ al are 
the absolute monochromatic fluxes for the standard star, 
then 

fx = qx x /-ii 0 -o.4 (m — mcaI ) ) (3 ) 

where m is the observed magnitude in a given band for 
which A is, for example, the mean wavelength of the filter. 
The q\ (~ 1) factors are required since the filters have 
finite widths and are not usually narrow in the infrared. 
It would be ideal to calculate the q\ factors from spec- 
trophotometric observations but in practice models are 
adopted (see the appendix in AAM94a for a discussion 
on this particular point). 

For the IRFM, the absolute flux calibration for the 
standard star must be done in such a way that the tem¬ 
peratures obtained in different wavelengths are close to 
each other and, when properly averaged, consistent with 
the direct temperatures. A calibration of this kind is that 
of AAM94a and is the one we have used in our work. 
Even though it was derived using direct temperatures 
of giant stars, this calibration provides good tempera¬ 
tures for dwarf stars as well, as it was shown in RM04b, 
where the comparison of IRFM angular diameters with 
measured ones for 7 main sequence stars showed good 
agreement (see also EH for an update). 

The models we have adopted for the ipx fluxes are Ku- 
rucz SEDs (Spectral Energy Distributions) as given in 
the library of theoretical spectra by Lejeune et al. (1997). 
From them and the response functions of the J, H, K 
and L filters of the TCS system (AAM94b, AAM98) we 
calculated grids of q\ and A t heo factors in the ranges 
3500 K< T eff < 9000 K, -3.5 < [Fe/H] < +0.5, and 
0.0 < logg < 5.0, where Kurucz models are available. 
Figure [I] illustrates the dependence of the temperature 
indicator in the IRFM (the q\ x i? theo factors, see (0J as 
a function of atmospheric parameters and filter. 

In the ranges of most interest to us (T e ff > 3800 K for 
giants, T e ff > 4500 K for dwarfs), the use of more sophis¬ 
ticated models will not have a substantial effect on the 
IRFM given that the model dependence of the method is 
relatively small. Also, unlike the UV and optical regions 
of the spectrum, the density of lines in the infrared is 
relatively low. The SEDs we took from Lejeune et al. 
(1997) are originally from Kurucz (1992). They were 
synthesized with Kurucz multi-million line list and at¬ 
mosphere models with scaled solar abundances and an 
approximate treatment of convective overshooting. The 
models were computed with a microturbulent velocity of 
2 km s -1 . 

The use of classical model atmospheres may not be a 
big issue for the IRFM. It is remarkable that the three- 
dimensional hydrodynamical simulations suggest that 
the effective temperature scale derived with the IRFM 
using one-dimensional hydrostatic models is the same as 
the one obtained with hydrodynamical models for turn¬ 


off stars with —3 < [Fe/H] < 0, and that for stars in 
the same metallicity range and approximately solar tem¬ 
perature the difference is only about 20 K (Asplund & 
Garcia Perez 2001). 

3. SAMPLE, PHOTOMETRY AND ATMOSPHERIC 
PARAMETERS ADOPTED 

Our sample consists of stars in the ranges of atmo¬ 
spheric parameters for which few data points were avail¬ 
able in AAM work, like the metal-poor and metal-rich 
extremes of the regions around T e ff ~ 5800 K and be¬ 
low 5000 K. The planet-hosting star sample is a metal- 
rich one given the well stated metallicity-enhancement of 
stars with planets (see e.g. Gonzalez et al. 2001, SIM04). 
Other metal-rich stars were taken from Castro et al. 
(1997) and Pompeia et al. (2002). The very metal-poor 
stars added are mainly from Ryan et al. (1991, 1996), 
Norris et al. (2001, mostly hot FG unevolved stars), 
and McWilliam et al. (1995); while the metal-poor cool 
dwarfs are from Yong & Lambert (2003a, 2003b). 

The distribution of our sample stars in the T e s vs. 
[Fe/H] plane is given in Fig. Q 

The distribution of the combined sample (i.e. including 
AAM stars) in T e ff and [Fe/H] is shown in Fig. [3] The to¬ 
tal sample is shown with the solid histogram, the dotted 
histogram is the AAM sample. Some of the peaks in the 
giant [Fe/H] distribution correspond to giants in clusters 
included in AAM work while those at [Fe/H] = 0.0, —1.0, 
and —1.5 in the dwarf [Fe/H] distribution are the stars 
with kinematic metallicity assignment (see 43.211 . Note 
how the dwarf metal-rich extreme is now much more pop¬ 
ulated, mainly with the planet-hosting stars. The very 
metal-poor end of both dwarfs and giants is also bet¬ 
ter represented, although still in small numbers. Some 
hot stars (T e ff > 7500 K) are also included but they are 
almost in the limit of applicability of the IRFM. 

3.1. Photometry and reddening corrections 

The infrared photometry for the AAM sample was 
taken from AAM94b and AAM98. This photometry is 
given in the TCS system and has a mean accuracy of 0.02 
mag. For our sample, Johnson and 2MASS photometry 
was adopted and transformed into the TCS system by 
means of the transformation equations of AAM94b and 
AAM98 (Johnson to TCS) and the following transfor¬ 
mation equations (2MASS to TCS) we derived for both 
dwarfs: 

•Arcs = •Amass + 0.001 - 0.049( J - A') 2M ass , (4) 
Atcs = AAmass — 0.018 + 0.003(J — Admass , (5) 
A'tcs = Admass - 0.014 + 0.034( J - A') 2M ass , (6) 

with standard deviations of a = 0.038 (N = 104), er = 
0.039 (N = 103) and a = 0.035 (TV = 107), respectively 
(N is the number of stars in every fit); and giants: 

•Arcs = A 2 MASS + 0.012 - 0.063( J — A') 2M ass , (7) 
AAres = AA 2 mass - 0.013 - 0.009(J — A') 2M ass , (8) 
A'tcs = Admass - 0.014 + 0.027(J - A') 2M ass , (9) 

with standard deviations of er = 0.036 (N = 43), cr = 
0.038 (N = 40) and a = 0.027 (N = 43), respectively. 

The V magnitudes were mostly obtained from the dif¬ 
ferent catalogs included in the General Catalogue of Pho¬ 
tometric Data (Mermilliod et al. 1997, hereafter GCPD). 
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Fig. 1.— Theoretical qR factors adopted as a function of effective temperature for typical atmospheric parameters of dwarf flog g = 4.5), 
giant (logg = 2.5), metal-rich ([Fe/H] = 0.0), and metal-poor ([Fe/H] = —2.0) stars; and filter. The units of Atheo are A for log g = 2.5. 
For log g = 4.5 the lines have been shifted by —2.0 dex. A rough estimate leads to qR x 7/)'/, almost independently of the band. 



Fig. 2.— Distribution of our sample stars in the T e g vs. [Fe/H] 
plane, open circles: dwarf stars; filled circles: giant stars. 


For the globular and open clusters we adopted the red¬ 
dening corrections given by Kraft & Ivans (2003) and 
AAM99a, respectively. For the field stars in AAM sam¬ 
ple, their E{B — V) values were adopted, while for our 
sample stars the reddening was estimated as explained 
below. 

We estimated E(B — V) values from the maps of 
Sclilegel et al. (1998) and Burstein & Heiles (1982), sev¬ 
eral extinction surveys included in Hakkila et al. (1997) 
code, and empirical laws by Bond (1980) and Chen et al. 
(1998). For stars closer than 75 pc, E{B — V) = 0 was 
adopted. The Arenou et al. (1992) extinction model in¬ 
cluded in Hakkila et al. code seems to systematically 
overestimate the reddening of stars with d < 0.5 kpc, 


TABLE 1 

Extinction Ratios Adopted, k = E(co\or)/E(B — V) or 
A hand /E(B - V). 


Color or band 

k (mag) 

(V - A'tcs) 

2.74 

(V - Kj) 

2.72 

Arcs 

0.84 

Arcs 

0.56 

A'tcs 

0.36 

Kj 

0.38 

Atcs 

0.15 


and so we have given a lower weight to the E{B — V ) 
values obtained from Arenou et al. map for stars closer 
than 500 pc. A latitude-dependent weight was adopted 
for the reddening obtained from the Schlegel et al. and 
Burstein & Heiles maps, in such a way that the high lat¬ 
itude objects (\l\ > 50) had a larger weight and the low 
latitude objects (\l\ < 30) a weight that was close to zero. 

There is evidence that the E(B — V) values obtained 
from the Schlegel et al. map, which is based on dust maps 
by DIRBE/COBE observations, overestimates the red¬ 
dening, so we have multiplied the Schlegel et al. E(B—V) 
values by 0.8 (Arce & Goodman 1999; Chen et al. 1999; 
Beers et al. 2002; Dutra et al. 2002, 2003). 

The extinction ratios k = E (color)/ E(B — V), or k = 
Aband /E(B — V ) given in Table |T| were obtained from 
Schlegel et al. (1998), using the appropriate effective 
wavelengths of the filters. 

3.2. Metallicities 

Our basic reference for the [Fe/H] values is the Cayrel 
de Strobel et al. (2001) catalog. We have updated 
the catalog with recent spectroscopic determinations, 
the most important of them being Mishenina & Kov- 
tyukh (2001), Santos et al. (2001), Heiter & Luck 
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Fig. 3.— Distribution of the AAM (dotted lines) and combined (AAM and this work, solid lines) giant (top panels) and dwarf (bottom 
panels) samples in [Fe/H] and T e fj. N is the number of stars in 0.1 dex bins in [Fe/H] or 200 K bins in T e g. 


(2003), Stephens & Boesgaarcl (2002), Takada-Hiday 
et al. (2002), and Yong & Lambert (2003a). From now 
on, we will refer to this updated catalog as C03. 

For the stars with as many or more than four entries in 
the catalog, the standard deviation of the different mea¬ 
surements was adopted as the error bar for [Fe/H], which 
was assumed to be the average of the values reported in 
C03. For those with less than four entries a mean error 
of 0.15 dex was adopted. 

Based on the C03 catalog and Stromgren photometry, 
we were able to calculate the metallicities of ~ 20 % of 
the sample dwarf stars from the following photometric 
calibration: 

For 0.19 < (6 — y) < 0.35, with cr = 0.17 dex: 

[Fe/H] = —4.29 — 66.0mi + 444.2mi(6 — y) 

—782.4mi(6 — y) 2 (10) 

+(0.966 — 37.8mi — 1.707ci) log 77 , 

where 77 = mi — [0.40 — 3.0(6 — y) + 5.6(6 — y) 2 ]. 

For 0.35 < (6 — y) < 0.50, with cr = 0.13 dex: 

[Fe/H] = -3.864 + 48.6mi - 108.5m? 

—85.2mi(6 — y) + 190.6m?(6 — y) (11) 
+ [15.7mi - ll.lci + 17.7(6 - y)\ Cl . 

For 0.50 < (6 — y) 0 < 0.80, with a = 0.15 dex: 

[Fe/H] = -2.63 + 26.0mi - 41.3m? - 45.4mi(6 - y) 

+74.0m?(6 - y) + 17.0miCi . (12) 

A mean error of 0.2 dex was adopted for these pho¬ 
tometric metallicities. These formulae follow Schuster 
& Nissen (1989) parameterization, but our relations are 
based on a much larger number of stars (~ 800) and 


they extend to late K stars. The formulae are valid for 
log g > 3.4 and —2.5 < [Fe/H] < 0.4. For F dwarfs they 
extend to [Fe/H] ~ —3.5. The major part of the disper¬ 
sion comes from the metal-poor stars, for [Fe/H] > —1.5 
the standard deviations reduce to 0.1 dex. 

The metallicities derived by SIM04 were found consis¬ 
tent with C03 and were thus adopted for the sample of 
planet-hosting stars. Although the temperatures from 
the IRFM and those derived from a purely spectroscopic 
approach differ by 100 K the metallicities are not severely 
affected if another T e g scale is adopted (J. Melendez & 
I. Ramirez 2005, in preparation). 

The photometric metallicities derived by AAM96a and 
AAM99a were also adopted, when necessary, and a mean 
error of 0.25 dex was assumed. Kinematic metallicities 
were assigned as a last resource. For these stars (~ 10% 
of the dwarf and held giant samples) a mean error of 
0.5 dex in [Fe/H] was adopted. The criteria was [Fe/H] = 
0.0 for thin disk stars, [Fe/H] = —1.0 for thick disk stars, 
and [Fe/H] = —1.5 for halo stars. The temperatures for 
these stars are very unreliable and should be taken only 
as rough estimates. 

Most of the metallicites for the giant stars in clusters 
have been taken from Kraft & Ivans (2003). For the 
clusters not included in Kraft & Ivans work, AAM99a 
values were adopted (see AAM99a for the references). 

3.3. Surface gravities 

As we did with the metallicities, we adopted the mean 
values from the C03 catalog when as many or more than 
four entries were available. These include both spectro¬ 
scopic and Hipparcos gravities and even though it has 
been suggested that a systematic difference of the order 
of 0.1 dex exists between them, it has almost no effect 
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Fig. 4.— Comparison of AAM95 bolometric flux calibration 
with B98 theoretical BC scale and their effect on the Tiff scale. 
Upper panel: fluxes from the AAM95 calibration minus B98. The 
dotted lines show how the fluxes would need to be off as to produce 
an IRFM T e ff scale shifted by d=50 K. Lower panel: our IRFM 
temperatures minus B98 (see the text for details). 
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Fig. 5.— Difference between the fluxes from AAM99a calibration 
and the mean of the measurements found in the literature (filled 
circles) and M03 (open circles) as a function of T e ff and [Fe/H]. 
The dotted lines show how the fluxes would need to be off as to 
produce an IRFM T e ff scale shifted by d=50 K. 


on the T e ff derived with the IRFM (see e.g. Fig. 4 in 
AAM96a). For the stars with less than four entries a 
mean error of 0.30 dex was adopted. 

The exact value of log g is not required for the IRFM, a 
distinction between giant and dwarf is, in principle, rea¬ 
sonable. Thus, when no spectroscopic or Hipparcos grav¬ 
ities were available, the log g value was estimated from 
the spectral type and luminosity class as given in Sim- 
bad and the representative log g values given by Straizys 
(1995). An error bar of 0.6 dex was assumed in this case. 
The log g values of AAM96a and AAM99a were used for 
the remaining stars (0.7 dex of mean error). 

The surface gravities of the cluster giants were cal¬ 
culated from their distance moduli (Kraft & Ivans 
2003) and the bolometric corrections of AAM99a, tak¬ 
ing 0.85Mq as an estimate of the stellar masses. 

3.4. Bolometric fluxes 

Bolometric fluxes have been calculated from the K 
magnitudes and (K — K) colors, using the calibrations 
of AAM95 and AAM99a, which are based on integrated 
UBVRIJHK photometry. These calibrations have a 
mean accuracy of 1% and are in reasonable agreement 
with other flux determinations (AAM99a, RM04b). 

Systematically inaccurate fluxes lead to offsets and/or 
trends in the temperature scale. As an illustration, 
we have calculated self-consistent bolometric fluxes and 
temperatures from the theoretical bolometric correction 
(BC) scale of B98 and compared them to ours for 13 
dwarf stars 4 in Fig. 0 We derived /boi values from their 
tables of BC by calculating an initial flux /boi(l) from 
our IRFM T e ff and the adopted parameters. Then, us¬ 
ing the /boi from the calibration of AAM95 and /boi(l) 
we got T eff (2) = T e ff(/boi(l)//boi ) 1/4 and iterated. This 
procedure is safe, because it converges simultaneously to 
/boi and T e ff independently of the input values, which 
is not true if we use the measured angular diameters to 
iterate in T e g. 

FigureOlshows how different bolometric flux scales lead 
to different temperature scales. Lower temperatures are 
clearly correlated to lower fluxes, although the correla¬ 
tion is not linear. While the difference in the fluxes 

4 These 13 stars are used in rmto compare measured diameters 
and direct temperatures with the IRFM (see Tableful. 


increases by 7% (from ~ —5.5% at T e ff ~ 4000 K to 
~ 1.5% at T e ff ~ 6500 K), the difference in the tem¬ 
perature scales increases only by about 2% (from —60 

K or-1.5% at T eff ~ 4000 K to 20 K or ~ 0.5% at 

T e ff ~ 6500 K ). Also shown in Fig. Q] as dotted lines are 
the flux differences that would produce a T e g scale offset 
of ±50 K in the IRFM. 5 Note how this explains almost 
completely the differences found with B98. 

The comparison of AAM95 fluxes with B98 may be an 
extreme case (see RM04b) so even if the adopted bolo¬ 
metric fluxes are inaccurate, a conservative estimate is 
that the IRFM temperatures are reliable within 50 K of 
systematic error. In RM04b we showed that the mea¬ 
surements of bolometric fluxes are within ±1% of the 
fluxes from AAM95 calibration, so the typical errors due 
to the fluxes are probably around 20 K. 

The situation with the calibration for giants (AAM99a) 
is similar, as it is shown in Fig.|3 where the fluxes from 
the calibration for about 30 stars 6 are compared to the 
mean of the measurements we found in the literature 
(see a and the fluxes by Mozurkewich et al. (2003, 
hereafter M03). As before, we conclude that the effect of 
the flux calibration on the T e ff scale may not be larger 
than 50 K of systematic error, and that typically it should 
be around 25 K. 

4. IRFM TEMPERATURES AND DIAMETERS 
4.1. Determination ofT e g and 6 

From the definitions of the q\ and Rtheo factors it is 
shown that 

gA^theo = f cal 10 —L(m-m cal ) ’ 4 3 ) 

where the left hand side is calculated from the models 
and filter transmission functions while the right hand 
side is obtained entirely from observations (adopting a 
calibration for the bolometric fluxes). The qR factors 
are thus considered as the temperature indicators of the 
IRFM. 

5 This is only an approximation but it is useful for the present 
purpose, the two lines are: A/| )n ] = 1 — [(T e ff ± 50 K)/T e ff] 2,7 . 

6 These stars are used in o to compare measured diameters 
and direct temperatures with the IRFM (see Table HI. 
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(V-K) 

Fig. 6.— A(qR)% is the error in the qR factors due to mean 
photometric uncertainties of 0.02 mag in ( V — K ) and 0.03 mag in 
(m — K), where m = J,H,K,L (see Ea. I14> . A(qR)%/2.7 is an 
estimate of the corresponding error in T e ff. Solid line: [Fe/H] =0.0, 
dotted line: [Fe/H] = —1.0, dashed line: [Fe/H] = —2.0. The set of 
lines in the bottom are for giant stars. The lines at top and those 
that diverge are for dwarfs. 


In this work, the bolometric flux calibration adopted 
has the functional form /boi = 10~ OAK <&(V — K). Then, 
from Eq. urn we deduce the relative error in the ob¬ 
served qR factor: 


A {qR) 
{qR) 


[0.921A(ro - K)Y 



1/2 


(14) 


Note that A$ is a function of (V — K) and its error 
A(V — K). Then, adopting mean errors of A(V — K) = 
0.02 and A (m — K) = 0.03 we obtained the approximate 
uncertainties A {qR)% plotted in Fig. El 

Since, roughly, qR oc T 2 7 independently of the band 
(Fig. |TJ, the A(qR)%/2.7 given in Fig. El is a reasonable 
estimate of the uncertainty in T e ft introduced by photo¬ 
metric errors. Thus, in our implementation of the IRFM 
we have adopted a mean error of 3% in qR as the error 
due to the photometry. 

The errors introduced by the atmospheric parameters 
[Fe/H] and logg were explicitly considered since their 
updated values are a major improvement to AAM work. 

In the best cases four temperatures were derived, one 
for each of the bands: J, i7, K , and L. Only the first 
three of them were available in the dwarf sample. The 
sensitivity of {qR)j (the qR factors in the J band) to 
T e g decreases for T < 5000 K for both dwarf and giants 
(Fig. |TJ and thus no Tj (temperature derived from the 
J band) below 5000 K was used to get the final T e g. The 
same reason restricted the use of Th to Th > 4200 K for 
the dwarf sample. Also, not all the giant stars have L 
magnitudes measured and so not all of them have a Tl . 

Figured shows that, for our sample stars, the temper¬ 
atures derived in the different bands are consistent. The 
dispersion is larger for the dwarf sample but no trends 
are observed. The same result is obtained when the up¬ 
dated temperatures for the AAM stars are included. 

The adopted IRFM effective temperature and its error 


was obtained from: 

_ EJTKATQ- 1 ] N 

6ff Ei(A TO" 1 ’ eff E,(A TO" 1 ’ 

( 15 ) 

where i = J, H, K, L ; and N is the number of tempera¬ 
tures available. 

The angular diameters were calculated from Eq. CD 
using the adopted bolometric fluxes (from AAM95 or 
AAM99a calibrations) and the IRFM temperatures. 
Their errors were estimated by propagating a mean error 
of 1.5 % in the bolometric fluxes and the errors in T e g. 

Table 0 contains our sample stars, their adopted at¬ 
mospheric parameters logg, [Fe/H], fluxes from AAM95 
or AAM99a calibrations, reddening corrections, {V — K) 
colors in the TCS system (transformed from either John¬ 
son or 2MASS photometry), IRFM temperatures and 
their errors, and angular diameters and their errors. The 
table has been ordered alpha-numerically by the star 
standard names: HD numbers have been preferred, then 
BD and then G numbers. Stars with other nomencla¬ 
tures have also been ordered alpha-numerically. A more 
detailed version of this table as well as tables for the 
combined sample, including the updated T e g and atmo¬ 
spheric parameters of the AAM stars, are available upon 
request. 

4.2. The Alonso et al. stars 

By comparing the effective temperatures we derived for 
the AAM sample with the AAM96a and AAM99a tem¬ 
peratures (Fig. El we find good agreement (within error 
bars) and thus show that the present T e g scale is consis¬ 
tent with that by AAM within the ranges in common. A 
closer inspection of Fig. 0 reveals very small systematic 
differences of about 40 K (0.7% at solar temperature) for 
dwarf stars and 25 K (0.6% at 4500 K) for giants, which 
are, however, well within the uncertainties of our results 
(1.3%). The scatter is due to the random nature of the 
differences between our input data and that adopted by 
AAM, and the extreme cases can be easily identified by 
comparing the two sets of input data. Less than 9 % 
(20 %) of the giant (dwarf) stars in Fig. 0 exhibit differ¬ 
ences larger than our mean uncertainty of 60 K (75 K). 
If we take into account the uncertainties in AAM tem¬ 
peratures these numbers are even smaller. 

4.3. The need for improvement and extension 

One important improvement to the IRFM T e g scale 
is the use of more accurate input parameters. Even 
though individually the temperature differences between 
the temperatures we obtained and those by AAM are rel¬ 
atively small, overall they greately help to reduce the dis¬ 
persion in the T e g : color : [Fe/H] relations (Part II), thus 
allowing to separate real effects from random errors. The 
revised IRFM temperature scale may then have a better 
internal accuracy, which is a very useful improvement, 
particularly in the light of the need to reduce (or even 
better, recognize) the sources of uncertainty in chemi¬ 
cal abundance studies. With the internal uncertainty 
reduced, the systematics can be better quantified. 

Also important is the effort to include stars represen¬ 
tative of regions of particular interest in the atmospheric 
parameters space, like the metal-rich (e.g. for planetary- 
host star studies) and metal-poor (e.g. for globular clus¬ 
ter and halo star studies) extremes. As an example of 






Ramirez & Melendez 



Fig. 7.— Left panel: temperatures obtained from the J and L bands ( Tj , T /,) compared to those obtained from the K band (T/<-); open 
symbols: dwarfs, filled symbols: giants. Right panel: as in the left panel for Tj. j vs. T ^. 


TABLE 2 

Input Data Adopted, Effective Temperatures and Angular Diameters Derived in This Work. 


ID 

log g 

[Fe/H] 

/bol 

E(B-V) 

(V-K) 

T e g (K) 


6 (mas) 

cr(6») 

BD -09 1261 

4.00 

-0.93 

2.034E-09 

0.000 

2.923 

3967 

148 

0.157 

0.012 

BD -10 3166 

4.41 

0.45 

2.764E-09 

0.010 

1.953 

5228 

67 

0.105 

0.003 

BD -13 3442 

3.92 

-2.95 

2.381E-09 

0.000 

1.226 

6442 

80 

0.064 

0.002 

BPS CS 16085-0050 

1.00 

-3.10 

4.759E-10 

0.010 

2.225 

4876 

55 

0.050 

0.001 

BPS CS 22166-0030 

4.00 

-3.36 

1.239E-10 

0.033 

1.461 

6209 

78 

0.016 

0.000 


the importance of these selected regions we show in Fig. 
El the T e g vs. (V — K) relation defined by our sample 
stars as well as AAM96b relations for [Fe/H] = 0 and 
[Fe/H] = —3 (photometry sources and reddening correc¬ 
tions employed are discussed in detail in Part II). The 
AAM96b relation for metal-poor cool dwarfs seems to 
be in agreement with ours (even if the AAM96b relation 
is slightly extrapolated) but the corresponding AAM96b 
relation for hot metal-poor stars differs from ours. The 
‘redder’ points near 6500 K on the dwarf T e g vs. (V — K) 
plane are all very metal-poor stars ([Fe/H] ~ — 3). This 
is the place where our results show the most significant 
deviations with respect to AAM temperatures. In Part II 
we show that the synthetic T e g scales based on both 
MARCS and Kurucz models for metal-poor stars are con¬ 
sistent with our temperatures but not with AAM96b in 
this region. The inclusion of these stars allows to better 
define the trends of the IRFM temperature scale in that 
region, which has proven to be of fundamental impor¬ 
tance, for example, in the controversy of the primordial 
lithium abundance (Melendez & Ramirez 2004). 

5. THE IRFM COMPARED WITH OTHER METHODS 

Since we have shown that our stars effective tempera¬ 
tures are essentially consistent with the AAM T e g scale, 
in the comparisons below, the stars with IRFM tem¬ 
peratures are those for the whole sample (including the 
AAM stars), for which we have homogeneously calcu¬ 
lated IRFM temperatures. 


5.1. Angular diameter measurements 

The IRFM was explicitly developed to solve for both 
T e g and 6 simultaneously from the observed and theo¬ 
retical data. One of the fundamental tests of validity of 
the IRFM is thus the comparison with measured angular 
diameters. 

Richichi & Percheron (2002, hereafter CHARM) have 
collected most of the published stellar angular diameters 
to date, measured mainly with lunar occultation (LO) 
and long-baseline interferometry (LBI) techniques. A 
large number of giant stars constitute the major part 
of the catalog and only a small fraction of main sequence 
stars is included. 

We have compared the angular diameters that re¬ 
sult from our IRFM implementation with those given 
in CHARM in Fig. (see also Table 0- The stars 
in CHARM have been grouped according to the tech¬ 
nique employed to derive their angular diameters. First, 
a comparison is made with the indirect (IND) diameters 
of Cohen et al. (1999), then with LBI, and finally with 
LO diameters. 

The IND diameters are in excellent agreement with the 
IRFM for 6 < 8 mas, where stars in common are found. 
These are diameters derived from absolutely calibrated 
stellar spectra in the infrared, which is the main reason 
for the good agreement. These are not direct measure¬ 
ments and will not be used to derive direct temperatures 
in the comparisons below. 
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Fig. 8. — Difference between the temperatures derived in this work (TW) and those given by AAM96a and AAM99a as a function of 
our T e ff values and adopted [Fe/H] for both dwarf (open circles) and giant (filled circles) stars. The dotted lines are the mean errors in our 
T efS values (75 K for dwarf and 60 K for giant stars.). 




Fig. 9.— T e ff vs. (V — K) relation for dwarfs as given by our 
sample stars (filled circles: [Fe/H] < —1.5; open circles: [Fe/H] > 
—1.5) and according to AAM96b temperature scale for [Fe/H] = 0 
(solid line) and [Fe/H] = —3 (dotted line). 


The LBI technique uses the interference of the wave- 
fronts from the starlight entering the telescope to resolve 
smaller angular sizes. For the stars with large diameters, 
more than one LBI entry is often found in CHARM. The 
mean value of those entries and its standard deviation 
for each star were adopted for the comparison (hence the 
large error bars of some points in Fig. [Til. otherwise the 
single value and its internal error were used. Comparison 
with the IRFM for 29 stars shows agreement within 1.5 
mas for 9 > 10 mas and better than 0.5 mas below that 
value. 

Less measurements with LO are available due to the 
limited sky coverage of the Moon (~ 10%) and other ob¬ 
servational difficulties. Nevertheless, the LO diameters 
are reliable and have been shown to be well scaled to the 
LBI ones. When compared to the IRFM, a very good 
agreement is also found. The differences with the five 
LO measurements (average in some cases) are less than 
0.5 mas in the range 3 < 9 < 20 mas. 

From the CHARM angular diameters and the bolomet- 
ric flux measurements of Blackwell et al. (1979, 1990), 
Blackwell & Lynas-Gray (1994, 1998), Leggett et al. 
(1986), and M03; we obtained the direct temperatures 
of the stars in common with our held giant sample and 
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TABLE 3 

Giant Stars With Measured Angular Diameters in CHARM and M03 



HD 

[Pe/H] 

nCHARM 

"LD 

1 

fl M03 

e LD 

rpIRFM 

J eff 

fhol 

A>ol^ 

A>ol 

0 IRFM ( mas ) 

T-.dir 5 
i eff 

(M03) 6 


3546 

-0.72 

1.77 

± 

0.08 

LBI 




4935 

± 

52 

5.81E-07 



1.71 

± 

0.05 

4857 

± 

111 



3627 

0.17 

4.21 

± 

0.05 

LBI 

4.14 

± 

0.04 

4343 

± 

45 

2.07E-06 

2.12E-06 

2.12E-06 

4.18 

± 

0.11 

4355 

± 

30 

4392 ± 

27 

3712 

-0.10 

5.58 

± 

0.13 

LBI 

5.61 

± 

0.06 

4553 

± 

48 

4.68E-06 

4.70E-06 

4.70E-06 

5.72 

± 

0.15 

4614 

± 

55 

4602 ± 

29 

6860 

-0.07 

14.08 

± 

0.76 

LBI 

13.75 

± 

0.14 

3824 

± 

41 

1.20E-05 

1.26E-05 

1.26E-05 

12.98 

± 

0.34 

3719 

± 

102 

3763 ± 

23 

9927 

-0.01 

3.76 

± 

0.07 

LBI 




4380 

± 

48 

1.56E-06 

1.60E-06 


3.57 

± 

0.09 

4292 

± 

43 



10380 

-0.27 

2.81 

± 

0.03 

LBI 




4132 

± 

46 

8.56E-07 



2.97 

± 

0.08 

4247 

± 

28 



12533 

-0.07 

7.73 

± 

0.20 

LBI 

7.81 

± 

0.08 

4259 

± 

45 

6.12E-06 

6.66E-06 

6.66E-06 

7.47 

± 

0.19 

4277 

± 

57 

4254 ± 

27 

12929 

-0.25 

6.87 

± 

0.04 

LBI 

6.83 

± 

0.07 

4501 

± 

50 

6.19E-06 

6.33E-06 

6.33E-06 

6.73 

± 

0.18 

4478 

± 

21 

4493 ± 

28 

18884 

0.00 





13.24 

± 

0.26 

3718 

± 

46 

9.42E-06 


9.57E-06 

12.17 

± 

0.35 




3578 ± 

37 

29139 

-0.18 

20.87 

± 

0.57 

LBI 

21.10 

± 

0.21 

3883 

± 

44 

3.06E-05 

3.34E-05 

3.33E-05 

20.10 

± 

0.55 

3895 

± 

56 

3871 ± 

24 

29139 

-0.18 

20.26 

± 

1.63 

LO 




3883 

± 

44 

3.06E-05 

3.35E-05 


20.10 

± 

0.55 

3956 

± 

160 



44478 

-0.05 

13.97 

± 

0.16 

LBI 

15.12 

± 

0.15 

3667 

± 

44 

1.11E-05 

1.12E-05 

1.12E-05 

13.55 

± 

0.38 

3623 

± 

25 

3483 ± 

22 

44478 

-0.05 

13.45 

± 

0.21 

LO 




3667 

± 

44 

1.11E-05 



13.55 

± 

0.38 

3684 

± 

32 



61338 

0.00 

3.07 

± 

0.78 

LO 




3869 

± 

48 

7.44E-07 



3.16 

± 

0.09 

3925 

± 

498 



62509 

-0.02 

7.98 

± 

0.06 

LBI 

7.98 

± 

0.08 

4833 

± 

50 

1.16E-05 

1.15E-05 

1.18E-05 

7.99 

± 

0.20 

4823 

± 

25 

4858 ± 

30 

62721 

-0.27 

3.01 

± 

0.00 

LBI 




3988 

± 

48 

6.94E-07 



2.87 

± 

0.08 

3894 

± 

15 



76294 

-0.01 

3.24 

± 

0.08 

LBI 




4817 

± 

50 

1.86E-06 



3.22 

± 

0.08 

4806 

± 

61 



80493 

-0.26 

7.69 

± 

0.26 

LBI 

7.54 

± 

0.08 

3851 

± 

42 

4.03E-06 

4.10E-06 

4.10E-06 

7.42 

± 

0.20 

3798 

± 

65 

3836 ± 

24 

94264 

-0.20 

2.54 

± 

0.03 

LBI 




4670 

± 

51 

1.05E-06 



2.58 

± 

0.07 

4701 

± 

33 



99998 

-0.39 

3.21 

± 

0.03 

LBI 




3919 

± 

45 

7.98E-07 



3.19 

± 

0.09 

3905 

± 

23 



99998 

-0.39 

3.20 

± 

0.51 

LO 




3919 

± 

45 

7.98E-07 



3.19 

± 

0.09 

3911 

± 

312 



102224 

-0.44 

3.23 

± 

0.02 

LBI 




4378 

± 

46 

1.36E-06 



3.33 

± 

0.09 

4447 

± 

22 



112142 

0.00 

5.90 

± 

0.70 

LO 




3647 

± 

47 

1.75E-06 



5.45 

± 

0.16 

3504 

± 

208 



113226 

0.11 

3.23 

± 

0.06 

LBI 

3.28 

± 

0.03 

5049 

± 

59 

2.25E-06 

2.23E-06 

2.21E-06 

3.23 

± 

0.09 

5035 

± 

47 

4981 ± 

31 

124897 

-0.55 

20.87 

± 

0.52 

LBI 

21.37 

± 

0.25 

4231 

± 

49 

4.69E-05 

4.89E-05 

4.86E-05 

20.95 

± 

0.58 

4284 

± 

56 

4226 ± 

29 

135722 

-0.40 

2.75 

± 

0.01 

LBI 

2.76 

± 

0.03 

4834 

± 

50 

1.35E-06 

1.37E-06 

1.41E-06 

2.72 

± 

0.07 

4831 

± 

20 

4851 ± 

32 

150997 

-0.28 

2.52 

± 

0.11 

LBI 

2.62 

± 

0.03 

4948 

± 

54 

1.26E-06 

1.27E-06 

1.26E-06 

2.52 

± 

0.07 

4947 

± 

111 

4841 ± 

36 

164058 

-0.15 

10.04 

± 

0.21 

LBI 

9.86 

± 

0.13 

3927 

± 

42 

8.29E-06 

8.38E-06 

8.40E-06 

10.23 

± 

0.27 

3974 

± 

45 

4013 ± 

30 

189319 

0.00 

6.09 

± 

1.35 

LBI 

6.23 

± 

0.06 

3877 

± 

41 

2.75E-06 

2.86E-06 

2.86E-06 

6.04 

± 

0.16 

3899 

± 

433 

3858 ± 

24 

197989 

-0.12 

5.23 

± 

0.60 

LBI 

4.61 

± 

0.05 

4710 

± 

52 

3.57E-06 

3.61E-06 

3.63E-06 

4.67 

± 

0.12 




4757 ± 

30 

198149 

-0.18 

2.65 

± 

0.04 

LBI 




4907 

± 

54 

1.38E-06 

1.42E-06 


2.67 

± 

0.07 

4960 

± 

42 



210745 

0.25 

5.31 

± 

0.11 

LBI 

5.23 

± 

0.05 

4482 

± 

51 

3.21E-06 

3.27E-06 

3.27E-06 

4.89 

± 

0.13 




4351 ± 

27 

214868 

-0.25 

2.63 

± 

0.05 

LBI 




4303 

± 

47 

7.69E-07 

7.59E-07 


2.60 

± 

0.07 

4260 

± 

44 



217906 

-0.11 

17.02 

± 

1.42 

LBI 

17.98 

± 

0.02 

3648 

± 

43 

1.50E-05 

1.55E-05 

1.52E-05 

15.93 

± 

0.45 

3557 

± 

149 

3448 ± 

13 

221115 

0.04 

1.61 

± 

0.17 

LBI 




4980 

± 

63 

4.79E-07 

4.89E-07 


1.53 

± 

0.05 

4877 

± 

258 



222107 

-0.50 

2.66 

± 

0.08 

LBI 




4605 

± 

49 

1.09E-06 



2.70 

± 

0.07 

4637 

± 

72 




Technique employed to derive ^L£>: LBI. long-baseline interferometry, LO. lunar occultation; as given in CHARM. 

^ Bolometric fluxes from the AAM99a calibration. Units are erg cm - ^ s — 

^ Mean of bolometric flux measurements in erg cm - ^ s — See the text for references. 

4 _Q _1 

Bolometric fluxes measured by M03 in erg cm s 

Direct temperatures from the CHARM angular diameters and the mean of the bolometric flux measurements (when available) or the flux from AAM95 calibration. A 
mean error of 1.5% was adopted for the fluxes. 

^ Direct temperatures from the M03 angular diameters and bolometric flux measurements. A mean error of 1.5% was adopted for the fluxes. 


compared them with those from the present IRFM im¬ 
plementation in Fig. El (see also TableEJ). When no flux 
measurements were available, the fluxes from AAM99a 
calibration were adopted. Although the error bars are 
larger than 200 K for some stars, we may safely state that 
a reasonable agreement is found in the ranges covered, 
i.e. 3500 K< T eS < 5000 IC and -0.7 < [Fe/H] < 0.2. 
The mean difference (T/g™ — T^g r ), with the direct tem¬ 
peratures from LBI measurements only, is 9 ± 58 K. The 
IRFM temperatures start to become larger than the di¬ 
rect ones at T e g < 3800 K. Discarding the points below 
that limit, the difference reduces to only 4 ± 57 K. The 
zero point of the IRFM T e g scale for giants may thus be 
the absolute one, but the reader should be cautioned that 
this result is currently only proved for stars in the ranges: 
3800 K< T eS < 5000 K and -0.7 < [Fe/H] < 0.2, where 
the comparison is made and no trends are found. 

A more recent interferometric work is that of M03, who 
determined the angular diameters of 85 stars, 19 of which 
are included in our field giant sample. These diameters 
allow for a better comparison than the previous one with 
CHARM given the homogeneity of M03 results. The 
comparison between our IRFM diameters and the mea¬ 
surements is shown in Fig. El (see also Table 0. Below 
9 ~ 10 mas the agreement is excellent but above that 
limit the IRFM diameters appear to be lower than those 
by M03. The reason for this is the very low tempera¬ 
tures of the latter group of stars. As we stated above, 
the IRFM produces high temperatures compared to the 
direct ones below 3800 K. Thus, we have plotted the stars 


with T e ff < 3800 K with open symbols in Fig. n2to show 
that they are the ones that depart most from the 1:1 line. 
Note also that M03 work uses optical interferometry and 
so their corrections from uniform-disk to limb-darkened 
diameters may be inaccurate due to uncertainties in the 
modeling of the TiO lines (Dyck & Nordgren 2002). 

M03 also measured bolometric fluxes and provide di¬ 
rect temperatures. They are compared to the IRFM tem¬ 
peratures of the 19 stars in common in Fig. El As we 
did before, the stars with T e g < 3800 K have been plot¬ 
ted with open symbols. It is more clear now that the 
IRFM temperatures exceed the direct ones for the very 
late-type giants. Excluding them, however, the agree¬ 
ment is excellent, as the mean difference (T/g™ — T^g) 
is only 8 ± 61 K. 

The excellent agreement found with the measured 
angular diameters and direct temperatures of the gi¬ 
ant stars is not completely surprising, given that the 
monochromatic flux calibration in the infrared was de¬ 
rived in such a way that the direct temperatures of a 
given number of giant stars were well reproduced by the 
IRFM (AAM94a). 

The situation for the main sequence stars is more com¬ 
plicated since the angular diameters of dwarf stars are 
very difficult to measure. Kervella et al. (2004) have col¬ 
lected uniform-disk angular diameters of several dwarfs 
and subgiants, all measured by interferometry, and have 
uniformly corrected them for limb darkening. 

In Fig. Elwe compare the IRFM diameters with those 
given by Kervella et al. (2004) for 10 dwarf stars and 
2 subgiants (see also Table EJ. The two subgiants are 
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Fig. 10.— Comparison between the angular diameters from 
CHARM (#ld) and the diameters from the IRFM (#irfm) for 
giant stars. Open circles: indirect (IND) diameters from Cohen 
et al. (1999) for 55 stars; filled circles: diameters from long base¬ 
line interferometry (LBI) for 29 stars; stars: lunar occultation (LO) 
measurements for 5 stars. Error bars are individual uncertainties 
or standard deviations, the latter when more than one measure¬ 
ment was available. The differences between the IRFM angular 
diameters and the IND diameters, LBI, and LO measurements are 
also shown with error bars containing the uncertainties in both $ld 
and #irfm- 


HD 121370 and HD 198149. The original sources for the 
measured angular diameters include Lane et al. (2001), 
Kervella et al. (2003), Nordgren et al. (1999, 2001), 
and M03; as well as some other papers in preparation 
(see Kervella et al. 2004). For Procyon (HD 61421) and 
HD 121370 more than one value is given by Kervella et al. 
(2004); only the more accurate VLTI measurement was 
considered. Also shown in Fig. rm and Table El is the 
comparison of the IRFM results with the transit obser¬ 
vations of Brown et al. (2001) for the planet-hosting star 
HD 209458. 

The comparison shows that the IRFM diameters are 
in excellent agreement with the measurements, to bet¬ 
ter than 0.2 mas in 9. Also, this result seems to be 
valid for the wide temperature range from 4000 K to 
6500 K, although it is restricted to almost solar metal- 
licity (Fig. 1151) . The stars in TableElspan the metallicity 
range from —0.54 to 0.25 dex. 

Table E| also contains the bolometric fluxes from 
AAM95 calibration and those measured by several au¬ 
thors. When a measured /boi was available, the direct 
temperature was derived from it and the measured an¬ 
gular diameter using Eq. 09 , otherwise the flux from 
AAM95 calibration was adopted. The mean of the mea- 
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Fig. 11. — Comparison between the IRFM temperatures and 
the direct ones (in K) for the giants in common with CHARM. 
Filled circles: LBI measurements, stars: LO measurements. The 
differences between the IRFM temperatures and the direct ones 
from CHARM are also plotted as a function of T e ff and [Fe/H] 
with error bars containing the uncertainties in both T e ff (IRFM) 
and T e ff (dir). The dashed line is at AT e ff = 4 K and the dotted 
lines at AT e ff = 4 zb 57 K. 
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Fig. 12.— Comparison between the angular diameters from M03 
($ld) an( f the diameters from the IRFM ($irfm) f° r 19 giant stars. 
Filled circles: T e ff > 3800 K, stars: T e ff < 3800 K. The differences 
between the IRFM angular diameters and those by M03 are also 
shown with error bars containing the uncertainties in both $ld an( f 
#IRFM • 
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TABLE 4 

Dwarf and Subgiant Stars With Measured Angular Diameters 


HD 

[Fe/H] 

$LD 0 

mas) 

T(K) 

/boi 1 

/boi 2 

References 2 

0IRFM ( 

'mas) 

(K) 

10700 

-0.54 

2.08 

± 

0.03 

5372 

± 

65 

1.163E-06 

1.154E-06 

1,2,3 

2.047 

± 

0.059 

5319 ± 43 

16160 

-0.03 

0.94 

± 

0.07 

4714 

± 

67 

1.687E-07 

1.725E-07 

2,3 

1.013 

± 

0.032 

4930 ± 185 

22049 

-0.12 

2.15 

± 

0.03 

5015 

± 

56 

1.024E-06 



2.204 

± 

0.059 

5078 ± 40 

26965 

-0.28 

1.65 

± 

0.06 

5068 

± 

63 

5.313E-07 

5.313E-07 

1,2,3 

1.555 

± 

0.045 

4920 ± 91 

61421 

-0.01 

5.48 

± 

0.05 

6591 

± 

73 

1.837E-05 

1.822E-05 

3,4,5,6,7,8 

5.405 

± 

0.145 

6532 ± 39 

88230 

-0.12 

1.29 

± 

0.04 

3950 

± 

161 

1.373E-07 



1.301 

± 

0.108 

3967 ± 63 

121370 

0.25 

2.20 

± 

0.03 

6038 

± 

75 

2.210E-06 

2.206E-06 

1,8 

2.234 

± 

0.065 

6081 ± 47 

128620 

0.20 

8.51 

± 

0.02 

5759 

± 

70 

2.677E-05 



8.548 

± 

0.244 

5771 ± 23 

128621 

0.20 

6.01 

± 

0.03 

5201 

± 

65 

8.653E-06 



5.957 

± 

0.173 

5178 ± 23 

131977 

0.09 

1.23 

± 

0.03 

4571 

± 

52 

1.921E-07 

2.010E-07 

2 

1.149 

± 

0.031 

4469 ± 57 

198149 

-0.18 

2.67 

± 

0.04 

4907 

± 

54 

1.377E-06 

1.414E-06 

1 

2.670 

± 

0.071 

4939 ± 41 

209100 

-0.02 

1.89 

± 

0.02 

4642 

± 

54 

5.069E-07 

5.000E-07 

2 

1.810 

± 

0.050 

4527 ± 29 

209458 

-0.01 

0.23 

± 

0.02 

5993 

± 

71 

2.279E-08 



0.230 

± 

0.006 

6049 ± 202 


j_ _n _1 

Bolometric fluxes from the AAM95 calibration. Units are erg cm s 
2 _o _ i 

Mean of bolometric flux measurements in erg cm s . References: 1. Blackwell &; Lynas-Gray (1994), 2. Blackwell &: Lynas-Gray (1998), 3. AAM95, 4. Code et al. 
(1976), 5. Beeckmans (1977), 6. Blackwell & Shallis (1977), 7. Smalley &; Dworetsky (1995), 8. M03. 
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Fig. 13.— Comparison between the IRFM temperatures and the 
direct ones (in K) for the giants in common with M03. Filled cir¬ 
cles: T e ff > 3800 K, stars: T e ff < 3800 K. The differences between 
the IRFM temperatures and the direct ones from M03 are also 
plotted as a function of T e ff and [Fe/H] with error bars containing 
the uncertainties in both T e ff (IRFM) and T e ff (dir). The dahed 
line is at AT e ff = 8 K and the dotted lines at AT e ff = 8 ± 61 K. 


surements of /boi is given in Table E| The direct temper¬ 
atures and their errors are also given in Tabled a mean 
error of 1.5% in /boi was assumed. 

The comparison between direct and IRFM tempera¬ 
tures is shown in Fig. The three points that depart 
most from the 1:1 line are HD 16160, HD 26965 and 
HD 209458. The uncertainty in the angular diameters 
measured for these stars is larger than 4% and, thus, 
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Fig. 14.— Comparison between measured angular diameters 
($ld) and diameters from the IRFM (#irfm) for the 13 dwarf (11) 
and subgiant (2) stars in Table The differences between the 
IRFM angular diameters and the measurements are also shown 
with error bars containing the uncertainties in both #ld and 
0IRFM- 
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Fig. 15.— Difference between the IRFM angular diameters and 
the measurements as a function of the metallicity and the IRFM 
effective temperature for the 13 dwarf (11) and subgiant (2) stars 
in Table |J] Units of A 6 are mas. The error bars contain the 
uncertainties in both $ld and #irfm- 
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Fig. 16.— Comparison between the IRFM temperatures and the 
direct ones (in K) for the 13 dwarf (11) and subgiant (2) stars in 
Table FTl The three open symbols are HD 16160, HD 26965, and 
HD 209458; they are not plotted in the residuals. The differences 
between the IRFM temperatures and the direct ones as a function 
of T e ff and [Fe/H] are also shown with error bars containing the 
uncertainties in both T e ff (IRFM) and T e ff (dir). The dahed line 
is at AT e ff = 18 K and the dotted lines at AT e ff = 18 =b 62 K. 


they are not useful to constraint any T e ff scale within 
200 K. Excluding these three points, the mean difference 
T m™_ T dir is 18±62 K. Consequently, the zero point of 
the IRFM T e ff scale for dwarfs is also in good agreement 
with the direct T e g scale, given that the 18 ± 62 K differ¬ 
ence is well within the average uncertainty of the IRFM 
temperatures (~ 75 K), at least for stars with temper¬ 
atures between 4000 K and 6500 K and metallicities in 
the range —0.55 < [Fe/H] < 0.25. 

5.2. The Sun, Procyon, and Arcturus 

The temperatures of these three stars are very well 
documented, although no general agreement concerning 
their exact values (the Sun being the exception) has been 
achieved, and should be used to test any T eS scale. 

Using the photometry given by AAM94a, the present 
IRFM temperature of the Sun results in 5756 ± 71 K 
which is 21 K lower than the standard value. AAM96a 
obtained 5763 ± 90 K and thus the difference with our 
solar T e ff may only be attributed to the details of the 
calculations. The fact that the IRFM does not exactly 
reproduce the solar temperature has been suggested to 
imply that the zero point of the IRFM T e g scale is off 
by some amount. This is by no means conclusive, as the 
photometry of the Sun is very uncertain and the errors in 
the observed magnitudes and colors are easily propagated 
to the temperature derived with this method. 

Procyon (HD 61421), whose spectral type, as given 


in SIMBAD, is F5IV-V, has an IRFM temperature of 
6591 ± 73 K. AAM96a obtained 6579 ± 100 K, the differ¬ 
ence is mainly due to a slightly different set of [Fe/H] 
and log g values adopted. Its direct temperature is 
6532 ±39 K, as given in Table 0J These different results 
may be consistent, considering the size of the error bars. 
Being a bright star, the photometry has a considerable 
uncertainty and may be the reason for the discrepancy, 
if any. Using the spectroscopic method, Takeda et al. 
(2002b, hereafter T02b) obtained 6576 K; Edvardsson 
et al. (1993, hereafter E93), on the other hand, ob¬ 
tained 6704 K from a theoretical calibration of the (b — y) 
color. A fit to spectrophotometric data by Castelli et al. 
(1997) resulted in T e g = 6650 K. Fitting of the Balmer 
line-profiles in the spectrum of Procyon leads to 6667 K, 
6470 K, and 6507 K according to Castelli et al. (1997), 
F98 and AP04, respectively. 

For the giant star Arcturus (HD 124897), the present 
IRFM implementation gives 4231 ± 49 K. That of 
AAM99a was 4233 ± 55 K. Its direct temperature, as 
given by Griffin & Lynas-Gray (1999) is 4290 ± 28 K; 
the value we derived for this star in Table [3] from several 
angular diameter and bolometric flux measurements is 
4284 ± 56 K, in perfect agreement with Griffin & Lynas- 
Gray result. M03 also derived a direct temperature for 
Arcturus, which is slightly lower (4226 K), but in better 
agreement with the IRFM. 

Although the IRFM temperature of Arcturus is lower 
by about 60 K, the temperature derived from several 
temperature-color calibrations based on the IRFM is 
closer to the direct one (4296 K, RM04a). Thus, the 
use of several of these calibrations allows to get rid of 
errors in single IRFM temperatures, which are mainly 
based on infrared magnitudes alone. We will come back 
to this point in Part II of this work. 

5.3. Other semi-direct methods 

In m we discussed some of the other methods devel¬ 
oped to derive T e g. As it is the case with the IRFM, 
they all introduce models to complement the available 
observations and so the differences that we will describe 
below are partly due to the different models adopted, al¬ 
though the observational input parameters may also be 
the reason of the discrepancies, if any. 

The comparison with the spectroscopic temperatures 
derived by three different groups is given in Fig. ^3 
In this and the next two figures the filled circles are 
dwarf stars, the open circles are held giants and the open 
squares are giants in clusters. 

T02b applied spectroscopic equilibrium to obtain the 
temperatures of 32 F-K dwarfs and subgiants, 23 of which 
are included in our dwarf sample. The comparison shows 
that the IRFM temperatures are higher by about 35 K 
but the dispersion is large (~ 100 K), mainly due to a 
metallicity effect. Below solar metallicity, T02b temper¬ 
atures are systematically lower; excluding the stars with 
[Fe/H] < 0, the mean difference reduces to 10 K, but 
then a clear trend with temperature is observed. 

Comparison with B03 spectroscopic temperatures of a 
sample of thin and thick disk stars suggests a mean dif¬ 
ference of —110 ±90 K, in which the dispersion is, again, 
mainly due to a trend of the differences with metal¬ 
licity. It seems that the differences reduce to zero at 
[Fe/H] ~ 0.25 and the B03 temperatures become sys- 
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Fig. 17. — Comparison of the IRFM temperatures with those 
derived by other authors. The figures show the IRFM temperatures 
minus those given by T02b, B03, and LT as a function of the IRFM 
T e ff and our adopted values of [Fe/H]. Filled circles: dwarfs; open 
circles: field giants; open squares: giants in clusters. 


tematically higher by 110 K as lower metallicities are 
reached. A similar difference was found in RM04b be¬ 
tween the IRFM temperatures and the spectroscopic 
ones of SIM04 for the planet-hosting stars but that one 
did not show any trend with metallicity, and only a slight 
trend with temperature, probably due to non-LTE effects 
and/or the bolometric fluxes adopted for the IRFM. 

In Fig. ^3 is also shown a comparison of the IRFM 
temperatures of 21 metal-poor field giants (—2.74 < 
[Fe/H] < —0.55) and 14 giants in the globular clusters 

M3 ([Fe/H]-1.5), M71 ([Fe/H]-0.8) and M92 

([Fe/H] 2.4) with those derived spectroscopically by 

Sneden et al. (1991), Kraft et al. (1992), and Shetrone 
(1996), which we refer jointly as LT (Lick-Texas group). 
The differences have a mean value of —65 ± 50 K, they 
show no large trends with T e g nor [Fe/H], and it is also 
independent of the population, i.e. it is basically the 
same for the giants in the field and in the three clusters. 
Thus, besides an offset of 65 K, the spectroscopic T e g 
scale of the metal-poor (—2.5 < [Fe/H] < —0.5) giants is 
differentially consistent with the IRFM. 

Using spectra synthesized from MARCS models, the 
filter transmission functions of the uvby-/3 system, and a 
set of stars to fix the zero point; E93 calculated a grid 
of synthetic colors from which the temperatures of their 
sample stars were determined. We found 50 stars (1 gi¬ 
ant) in common with our sample. The comparison with 
the IRFM is shown in Fig. where a mean difference 
of —75 ± 85 K is observed. No trend is apparent with 
metallicity, and probably a slight trend with temperature 
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Fig. 18.— As in Fig. 1171 for E93, F98, and AP04. 


is present, in the sense that the IRFM temperatures and 
those by E93 may be in agreement for T e ff > 6500 K but 
with a decreasing difference for lower temperatures. 

Two sets of temperatures based on the fitting of the 
Balmer line-profiles are shown in Fig.^1 Note how even 
though F98 and AP04 temperatures were derived using 
the same principle, the theoretical prescriptions adopted 
have an important effect on the final results. F98 tem¬ 
peratures are higher by about 100 K below 5800 K but 
slightly lower than the IRFM ones above 6000 K. There is 
apparent agreement at solar metallicity but disagreement 
both above and below [Fe/H] ~ 0. The solar tempera¬ 
ture is thus well reproduced by both the IRFM and F98. 
On the other hand, AP04 temperatures show very good 
agreement with ours in the wide temperature range from 
4800 K to 6600 K and it is likely that they are close to 
the IRFM ones even at low metallicites ([Fe/H] ~ —0.8). 

Note that the zero point of the E93 theoretical calibra¬ 
tion was set using a sample of stars with temperatures 
derived from Balmer line-profile fitting (Cayrel de Stro- 
bel 1990) and so the similarity between the differences 
found with E93 and F98 may be due to the same basic 
technique employed. We recall, however, that this tech¬ 
nique is very model dependent, as shown by the most 
recent and improved results of AP04, and that it is very 
likely that the results from the IRFM are in good agree¬ 
ment with those from the modeling of hydrogen lines. 

A method that was not discussed in m is the one 
that uses surface brightness vs. color relations based on 
measurements of stellar angular diameters (di Benedetto 
1998, hereafter dB98). Since the IRFM is also based on 
measured angular diameters, the zero points of both the 
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IRFM and the dB98 T e g scales should be the same. This 
is actually shown in Fig. CHI where the temperatures 
of 25 stars (12 giants) in common with our sample are 
compared to those given by dB98. The mean difference 
is — 1 ± 53 K, as expected, and the dispersion is fully 
consistent with the mean errors of the IRFM and dB98 
temperatures. The good agreement of the zero point is 
also true if dwarfs and giants are considered separately. 
Note, however, that the comparison with dB98 results is 
limited to [Fe/H] > —0.6, although the only point be¬ 
low —0.6 dex, a giant star with [Fe/H] = —1.64, shows a 
dB98 temperature consistent with the IRFM. 

Finally, the line-depth ratio technique has been applied 
to 181 dwarfs by K03, 44 of them are included in our sam¬ 
ple. It is claimed that very small errors, of the order of 
10 K or less, may be achieved with this technique if sev¬ 
eral calibrated ratios are used. The comparison with our 
IRFM temperatures, however, shows a dispersion that is 
larger than the one expected by the combination of the 
errors in the IRFM temperatures (~ 75 K) with those 
suggested by K03, which will be essentially equal to the 
uncertainty from the IRFM. It is true that very small in¬ 
dividual uncertainties (that is, for a single star) may be 
achieved with the line-depth ratio technique, but some 
systematic effects may enlarge the dispersion in the T e g 
scale derived with this technique, which is, however, po¬ 
tentially useful for studying intrinsic temperature fluc¬ 
tuations. K03 adjusted their T e g scale to reproduce 
the solar temperature, hence the good agreement with 
the IRFM around 5800 K and solar metallicity. From 
5500 K to 6200 K, however, K03 temperatures become 
lower than the IRFM ones. From 4500 K to 5500 K, a 
mean difference of about —90 K is observed, and it seems 
to go to zero again at lower temperatures. 

6. CONCLUSIONS 

Effective temperatures with a mean accuracy of about 
1.3% have been determined with the infrared flux method 
for a selected sample of stars as well as for the Alonso 
et al. (AAM) sample with the aim of calibrating the 


metallicity-dependent temperature vs. color relations in 
a companion publication. The status of this particular 
IRFM implementation is revised and it is shown that the 
IRFM temperatures are well scaled to the direct ones for 
main sequence stars. The case for giant stars was proved 
in the past and the present work confirms the good agree¬ 
ment. The present work improves the IRFM T e g scale 
by the use of up-to-date spectroscopic atmospheric pa¬ 
rameters and the inclusion of stars in the ranges where 
few or no data were previously studied. Our IRFM ap¬ 
proach, however, is largely based on the AAM study, 
and a comparison of our updated to AAM temperatures 
shows agreement within error bars. 

We point out that the influence of the bolometric flux 
calibrations adopted on our IRFM temperatures is im¬ 
portant because they may introduce systematic trends 
in the resulting T e g scale. However, we have also shown 
that the systematic trends may not be larger than 50 K, 
and so, if differences larger than 50 K are found with 
other T e g scales, they should not be attributed solely to 
the fluxes adopted in the IRFM. This number (50 K) is a 
conservative estimate, and thus even smaller differences 
may not be totally due to errors in the bolometric fluxes. 

Comparison of the stellar angular diameters from the 
IRFM with those measured by interferometry, lunar oc- 
cultation or transit observations shows good agreement. 
This work thus provides accurate stellar angular diam¬ 
eters, which are essential to calibrate extremely high- 
resolution long-baseline interferometric observations (e.g. 
VLT and CHARA interferometers). 

When combined with bolometric flux measurements, 
the interferometric angular diameters produce a temper¬ 
ature scale with a zero point that is only 10 K below 
that derived with the IRFM for the giant stars with 
[Fe/H] > —0.7 and temperatures above 3800 K and be¬ 
low 5000 K. The corresponding zero points of the dwarf 
T e g scales differ by only about 20 K (the IRFM being 
larger), and it is valid for [Fe/H] > —0.55 and temper¬ 
atures between 4000 K and 6500 K. No stars outside 
those ranges were available for comparison, and we can¬ 
not guarantee that the zero point of the IRFM T e g scale 
is essentially the absolute one everywhere in the atmo¬ 
spheric parameters space, as it is the case in the ranges 
mentioned above. 

The IRFM temperature of the Sun is only about 20 K 
below the standard value, and those of the well studied 
stars Procyon and Arcturus are consistent with the di¬ 
rect ones and those derived with other methods within 
the error bars. No perfect agreement is expected for these 
bright stars because the photometric uncertainties asso¬ 
ciated are large and easily propagated in the IRFM. Con¬ 
sequently, no empirical shifts to reproduce the tempera¬ 
ture of the Sun or other well studied stars were applied. 

The comparison with the temperatures obtained with 
other methods shows that the temperature scale of the 
lower main sequence is still unreliable at the 100 K level 
because larger discrepancies are found. Although con¬ 
sistent with the temperatures from Balmer line-profile 
fitting, the IRFM produces temperatures that are about 
100 K lower that those derived from the excitation equi¬ 
librium of Fe I, but exceptions are found. The temper¬ 
ature scale of giant stars is better defined, probably to 
a 50 K level or even better. The spectroscopic temper¬ 
atures of the metal-poor giants, both in the held and 
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in clusters, are about 65 K higher that those from the 
IRFM, but otherwise differentially consistent with them. 

With the IRFM temperatures described in this paper 
and the atmospheric parameters adopted, we will cali¬ 
brate a set of homogeneous temperature vs. color re¬ 
lations (Part II). The fact that the input parameters 
are somewhat inhomogeneous (several sources have been 
considered) will not have a significant effect on these re¬ 
lations but they may have some effect on the IRFM tem¬ 
peratures. 
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